Tunable Kondo-Luttinger systems far from equilibrium 
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We theoretically investigate the non-equilibrium current through a quantum dot coupled to one- 
dimensional electron leads, utilizing a controlled frequency- dependent renormalization group (RG) 
approach. We compute the non-equilibrium conductance for large bias voltages and address the 
interplay between decoherence, Kondo entanglement and Luttinger physics. The combined effect 
of large bias voltage and strong interactions in the leads, known to stabilize two-channel Kondo 
physics, leads to non-trivial modifications of the conductance. For weak interactions, we build an 
analogy to a dot coupled to helical edge states of two-dimensional topological insulators. 
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Understanding strongly correlated quantum systems 
far from equilibrium is an outstanding challenge in 
condensed-matter physics. Many of the theoretical ap- 
proaches that have been proven successful in treating 
strong correlations are inadequate once the system is 
driven out of equilibrium. Quantum dot devices pro- 
vide an ideal setting to study transport under non- 
equilibrium conditions, as they constitute comparatively 
simple model systems with high tunability |l|-|7j. Kondo 
physics plays a crucial role in understanding their trans- 
port properties H, H|. It has been shown that sev- 
eral effects in these devices will suppress or modify the 
Kondo screening, such as dissipation and the electron- 
electron interaction in Luttinger liquid quantum wires 
Io|-15j]. In this Letter, we study 




that couple to the dot 
non-equilibrium currents across quantum dots in nano- 
settings (see Fig. [^involving Kondo entanglement and 
Luttinger physics (16J. 

A quantum dot in the Kondo regime coupled to one- 
dimensional (Id) leads exhibits either a one-channel 
Kondo (1CK) or a two-channel Kondo (2CK) ground 
state dHHH, as the Luttinger parameter K is decreased; 
the control parameter corresponds to the interaction 
strength in the Id leads flij . The non-equilibrium prop- 
erties of this system were addressed only in an ex actly 
solvable limit [l7| or in the linear (low bias) region 18| . 
The full crossover in the non-equilibrium conductance 
between the 1CK and 2CK fixed points (l4j |. with much 
relevance to experiments, has not yet been addressed. In 
particular, interactions in Id wires are expected to result 
in a peculiar non-equilibrium transport 2^, 2l| . 

Here, we apply a non-equilibrium RG method to 
tackle these issues. We calculate the conductance for bias 
voltages large compared to the relevant Kondo scales. 
We identify signatures of intermediate 2CK behavior in 
the RG flow for all K < 1 which strongly modify the 
conductance profile. The low-temperature conductance 
is non-universal in the sense that it does not depend on 



FIG. 1: Setups considered in this Letter, (a) A quantum dot 
coupled to two Id Luttinger leads (L, R) where both left and 
right-moving electrons in the leads can be either in the spin 
up or spin down state, (b) A quantum dot coupled to (helical) 
edge states (1, 2) of topological insulators. The right and left 
moving electrons on the edges are spin polarized. 



V/T-zck only, where T2CK is the relevant 2CK scale. 

There is also a growing interest in Kondo physics in 
topological insulator (TI) systems 22h24| . Due to spin- 
orbit coupling, TIs have gapless helical edge states where 
the direction of the electron's spin and momentum are 
entangled. We shall extend our analysis to a quantum 
dot coupled to two helical edges of 2d TIs (see Fig. O^b)) 
where it has been shown that 2CK physics is stable even 
for weak repulsive electron-electron interactions [24| . 

Equilibrium properties. Let us focus on the setup of 
Fig. 1(a). We denote by gLR and Qll = 9rr the dimen- 
sionless inter-lead and intra-lead Kondo couplings, re- 



spectively n 
fixed points 




The RG analysis results in two infrared 
the 1CK and 2CK fixed points. 



In the former case, all Kondo couplings, gLR and gLL, 
are relevant under RG transformation and flows towards 
strong coupling, such that the two leads can be combined 
into a single effective lead. In contrast, the 2CK fixed 
point is reached when gLR remains small under RG, while 
9ll grows (and flow to intermediate coupling). Here, the 
two leads provide independent screening channels. This 
2CK fixed point is infrared stable for K < 1/2 (assuming 
that g LL 
For K- 
realized If 



9rr) 0113- 
= 1 (free electron leads), where 1CK physics is 
, 0] , the conductance reaches the unitary limit 



2 



at low temperatures, G(T) = 2G [1 -0((T/T K ) 2 )] where 
Go = e 2 /h is the conductance quantum. Here, Tk = 
Doe~ 1 ^ 9LIi+9LL '> is the Kondo temperature; Dq is an ul- 
traviolet cutoff (~ few Kelvins) whereas g LR and g\ L are 
the bare values of the coupling constants. For T 3> Tk, 
from G(T) oc g LR (T) 2 , one finds G(T) ~ 1/ In 2 (T/T K ). 

For all K < 1, gLR grows slower under RG than 
9ll(rr) ■ For K -C 1, one can solve the RG equa- 
tions analytically for large T. We may neglect gLR in 
the RG equation for gLL/RR to obtain the approximate 
solution gLL/RR,{T) l/ln(T/T£) with the shorthand 



K 



. The coupling gLR.(T) is found by sub- 



stituting the approximate solution for 9ll/rr{T) in the 
RG equation for the coupling <jfL_R. We evaluate: 




FIG. 2: (Color online) gLh(u) and gLn(u>) for various Lut- 
tinger parameters if. The bare couplings are g1 a , = g^ a = 
0.025, resulting in T* K ~ 4.2 x 10" 18 and ~ 2 x 10" 9 (in 
units of D = 1). The bias voltage is V = 0.6 > T K , T£. 
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where D* 



[g^ln^Do/T*)]^ 



with 1/s = (1/K - l)/2. 



We deduce that, for T » Tt, the conductance G(T) oc 
fl£ H (T) essentially follows T^ K - 1 / ln 4 (T/T£). Here, the 
power-law behavior is reminiscent of Luttinger physics 
whereas the logarithmic contribution is typical of Kondo 
correlations. Importantly, the conductance is not a uni- 
versal function of (T/T£) because transport arises from 
the sub-leading coupling gLR- 

For 1/2 < K < 1, the low-temperature physics is gov- 
erned by two scales, Tick < T 2 ck, with 1CK behavior 
for T <C Tick and 2CK behavior for Tick < T < T 2CK - 
In the limits K — > 1 and K — > 0, we have Tick = Tk and 
T 2 ck = T K , respectively, with Tk, T^ defined above. In 
general, Tick Tk due to interactions in the leads; also, 
T2CK > T K . In the presence of particle- hole symmetry, 
the conductance for 1/2 < K < 1 reaches the unitary 
limit as T —> 0. However, potential scattering is a rele- 
vant perturbation with a scaling dimension (l+K)/2 < 1 
and causes the conductance to decrease as T X I K ~ X as 
T — > 0; the leading irrelevant operator corresponds to 
the hopping (gLR) term between the two leads with scal- 
ing dimension (l/K + l)/2. Similarly, near the 2CK fixed 
point reached for K < 1/2 and g^L = gRR, the leading 
irrelevant operator (gLR term) has dimension 
and therefore one expects G(T) oc T X I K ~ 2 as T ->• 0. 

Non-equilibrium properties. We study the low- 
temperature conductance in the high-bias regime V 3> 
Tk,T K where the non-equilibrium RG method can be 
applied 0. (In the opposite limit, V <C T K , we expect 
the equilibrium results quoted above to be valid after 
replacing T — » V.) Setting ft = fee = e = 1, the non- 
equilibrium RG equations take the form: 
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where Q u = Q(D-\u+iT\) and /? = -1(+1) labels leads 
L(R). Further, T is the decoherence (dephasing) rate at 
finite bias which cuts off the RG flow Q : 



r = Tr^yW: (1 - [<wM] s 



(3) 



where the Fermi function obeys f a (oj) = 1/(1 + 
e (u)-na,)/Ty no ^ e that there exists an additional con- 
tribution to r from electron dephasing caused by a finite 
potential drop in the Luttinger liquid leads [20(, which 
will affect the subleading terms in T (given by gLR(V/2)). 
However, in the low-conductance regime of interest, this 
voltage drop is small and will be neglected henceforth. 
In general, the perturbative RG approach is valid for 
V > Tick, T 2CK - In the limit of V -> 0, Eqs. © reduce 
to the equilibrium RG equations (with the flow cut off 
by temperature), and we recover Eq. ([1]). 

The renormalized couplings are obtained by self- 
consistently solving Eqs. ([2]) and ((3]) [2]. As shown in 
Fig. [21 gLL{RR)(u) exhibit peaks for all values of K < 1, 
indicating that they grow under RG. For a given bias 
voltage, the Kondo coupling gLR.(u) shows a crossover 
from peak to dip structure as K decreases, traducing the 
fact that for a fixed bias voltage, 5la( w ) is either en- 
hanced or decreased compared to its bare value g® LR - Let 
us emphasize that for sufficiently large bias voltages, as 
soon as K < 1 the coupling gLR(u) exhibits a dip close 
to oj = ±V/2, signalling 2CK behavior. The singular be- 
havior at the peaks or dips is cut off by the decoherence 
rate (see Eq. (6)), while outside that regime the voltage 
serves to cut off the RG flow. 

From the Keldysh calculation up to second order in the 
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FIG. 3: (Color online) G(V) normalized to G = 3n(g° LR ) 2 /4: 
for various K and Kondo couplings as in Fig. 2. Here, "Equi- 
librium form" refers to a fit with the expression G(V) = 
3ngLR(<-o = 0) 2 /4 which consists in substituting qlr{i>j) by 
gLn(oJ = 0) for \u\ < V/2. For K = 0.5 we also show the an- 
alytical result from Eqs. (7,8) (stars). Inset: r(w = 0, V) for 
various K. The dashed lines are obtained from the analytical 
expression in Eq. (9). Here, Do = 1. 



FIG. 4: (Color online) G versus V/T2CK for K = 0.6 for 
various initial Kondo couplings g^ a i = gaa = 9 ( m units of 
Do = 1). 



For K < 1, we manage to obtain an approximate analyti- 
cal form for the couplings 9lr{u = V/2) and gLR^ = 0). 
Solving Eqs. ^ in the limit D — > 0, we find: 



tunneling amplitudes, the current reads: 



(L «-> R). (4) 



For small bare couplings g^ a , — <?° Q this perturbative 
calculation of / remains valid for V > Tj£, implying 
at high bias voltage V » T contributions to the cur- 
rent over a frequency window — V/2 < lu < V/2. For 
bias voltages T£ < V <C D , with decreasing V we find 
the differential conductance G(V) = dl/dV approaches 
the equilibrium form of the conductance G(T — > V) oc 
V 1 / K - 1 /\n\V/T^) with G(T) oc g LR (T) 2 (see Eq. Q 
and Fig. [3]). 

In the remainder, we analyze G(V) for larger bias volt- 
ages. For K = 1, we checked that the nonlinear conduc- 
tance satisfies G(V) oc l/ln 2 (^/T K ) for V > T K (H. 
Here, one can replace g L n{u) by g_Li?.(w = 0) w 5Lfl(r -> 
V). When decreasing if, the double peak structure in 
9lr{w) at w = ±V/2 turns progressively into dips which 
acquire a complex shape as a result of the decoherence 
rate Y and the electron-electron interaction which hin- 
ders the inter-lead electron tunneling. The effect be- 
comes more pronounced for small K values associated 
with the 2CK fixed point, rendering the "flat" approxi- 
mation <7l_r(oj) ~ 9lr(u = 0) not justified; see Fig. [3] 

To gain an analytical understanding of the small-if 
non-equilibrium regime, we may treat g R L{ui) within the 
interval -V/2 < lu < V/2 as a semi-ellipse [Tol ]. The 
current I reads 
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where gLR.(T — > V) is the equilibrium form of g^R in Eq. 
(P) with T replaced by V, and we have defined jF(K) = 
2 i+i(i-i/K)_ 1) with j:^ K = !) = x Using Eqs. (5) and 
(6), we obtain a closed expression for the conductance: 



K{V) 



(J) 



+ 12n[l--)W'(V), 



where W'(V) = dW/dV and: 



n{v) 

W(V) 



l/K 



vxrv/D* 2 )***- 1 ) 



(4 



(8) 



ln 4 [n//(T*) 



For completeness, we have kept the less dominant con- 
tribution - l/\n 5 (V/T£) in ft(V~). To rigorously define 
the function W(V), we need to provide an analytical ex- 
pression for the decoherence rate in Eq. (3). Using an 
analogous reasoning as for the non-equilibrium current I, 
to second order in g° LRl we extract 



4 In 4 



(9) 
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Close to K = 1/2, we can safely neglect contributions 
in {Qlr) 2 an d therefore to second order in g\ R , 

we find Y « (tt 2 / A)[g LR {uj — 0)] 2 . We have checked our 
analytical expression of Y against a numerical treatment 
of Eqs. (2) and (3); see inset in Fig. [31 Notably, the 
decoherence rate contributes to a "distinct" power law 



in the non-equilibrium conductance G(V) 



where {l/K 2 — l)/2 > 1/K — 1 for K < 1, rendering the 
second term in Eq. (7) to be subleading. The conduc- 
tance becomes smaller than its equilibrium counterpart 
since gLR,{w = V/2) < = 0). A comparison be- 

tween the analytical formula in Eq. (7) and the numeri- 
cal integration of Eqs. (2-4) is shown in Fig. [3] As our 
results are based on one-loop RG, we may expect both 
corrections to the power-law prefactors and further sub- 
leading terms upon including higher-loop contributions. 

Our results show that G(V) for voltages T 2 ck -C V <C 
Do is not an universal function of V/T2CK (even for fixed 
K): Fig. [4] displays G versus V/T2CK for various ini- 
tial Kondo couplings, with T2CK extracted from the RG 
flow. As it becomes also clear from Eqs. (7,8), the non- 
equilibrium conductance for V ^> T2CK is a function of 
both V/Dq and V/T£, and hence has a non- universal 
profile. This is again related to the fact that transport 
arises from the subleading coupling gLR- 

Topological insulators. We can extend these results to 
a quantum dot coupled to helical edges of 2d TIs; see Fig. 
1(b). In contrast with the setup of Fig. 1(a) where single- 
particle backscattering terms caused by the quantum dot 
cut the system into two separate parts, such backscatter- 
ing terms are now forbidden due to time-reversal symme- 
try of the helical edges [H, El- Following Refs. IM [27, 
the two edges of spinless electrons can be mapped onto 
spinful Luttingcr liquids with distinct Luttinger param- 
eters for charge (K c ) and spin (K s ) degrees of freedom: 
K c = K, and K s — 1/K with K being Luttinger param- 
eter for the helical edges. In equilibrium, the system is 
equivalent to an anisotropic 2CK model with the Kondo 



couplings g x J*p Z (a, (3 — 1,2) 24] and in contrast to the 



case in Fig. 1(a), the 2CK fixed point is stable for K < 1. 
In the limit of K — > 1~ where the anisotropy in the 
Kondo couplings is negligible, we have checked that the 
RG scaling equations for the Kondo couplings gn/22 and 
gi2 are identical to those of the Kondo couplings gLL/RR 
and g L R. We^ identify: gLL/RR -> 511/22, gLR -> 512, 
j7 — > K + 1/K — 1 [28] . Our results on non-equilibrium 
transport across a quantum dot coupled to weakly inter- 
acting Luttinger leads are directly applicable. 

Summary. We have studied non-equilibrium trans- 
port through a Kondo dot coupled to Luttinger-liquid 
leads and calculated the conductance profile at bias volt- 
ages larger than the Kondo scales of the system. The 
RG flow at large bias shows signatures of intermediate 
2CK physics for all Luttinger parameters K < 1. As 
the conductance G arises from the coupling gLR which 



is subleading, G(T — > 0) is not a universal function of 
V/T2CK as it also depends on V/D . Our results push 
forward the knowledge of correlation effects in nanosys- 
tems far from equilibrium and should stimulate further 
experimental works on transport through dots coupled 
to quantum wires and carbon nanotubes. We have also 
shown that our theoretical framework is applicable to the 
Kondo effect at the helical edges of topological insulators. 
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